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ABSTRACT 

We present simultaneous Chandra-High Energy Transmission Gratings (HETG) and Rossi X-ray Timing 
Explorer (RXTE) observations of a moderate flux 'soft state' of the black hole candidate 4U 1957+1 1. These 
spectra, having a minimally discernible hard X-ray excess, are an excellent test of modern disk atmosphere 
models that include the effects of black hole spin. The HETG data show, by modeling the broadband continuum 
and direct fitting of absorption edges, that the soft disk spectrum is only very mildly absorbed with Nh = 
1-2 X 10^' cm"^. These data additionally reveal AA13.449 Ne IX absorption consistent with the warm/hot 
phase of the interstellar medium. The fitted disk model implies a highly inclined disk around a low mass 
black hole rapidly rotating with normalized spin a* w 1. We show, however, that pure Schwarzschild black 
hole models describe the data extremely well, albeit with large disk atmosphere "color-correction" factors. 
Standard color-correction factors can be attained if one additionally incorporates mild Comptonization. We 
find that the Chandra observations do not uniquely determine spin, even with this otherwise extremely well- 
measured, nearly pure disk spectrum. Similarly, XMM-Newton/RXTE observations, taken only six weeks 
later, are equally unconstraining. This lack of constraint is partly driven by the unknown mass and unknown 
distance of 4U 1957+11; however, it is also driven by the limited bandpass of Chandra and XMM-Newton. 
We therefore present a series of 48 RXTE observations taken over the span of several years and at different 
brightness/hardness levels. These data prefer a spin of a* « 1, even when including a mild Comptonization 
component; however, they also show evolution of the disk atmosphere color-correction factors. If the rapid 
spin models with standard atmosphere color-correction factors of hj = 1 .7 are to be believed, then the RXTE 
observations predict that 4U 1957+1 1 can range from a 3 M0 black hole at lOkpc with a* « 0.83 to a 16 Mq 
black hole at 22 kpc with a* « 1, with the latter being statistically preferred at high formal significance. 
Subject headings: accretion, accretion disks - black hole physics - radiation mechanisms:thermal - X- 
rays:binaries 



1. INTRODUCTION 

Despite having an X-ray brightness comparable to and oc- 
casionally exceeding that of LMC X-1 or LMC X-3, and 
despite also being one of the few persistent black hole 
candidates (BHC), 4U 1957+11 has received relativ ely lit- 
tle at tention. Similar to LMC X-3 (see, e.g., .Wilms et al.l 
l2001h . 4U 1957+11 is usually in a s oft state that shows long 
term (h undreds of days) varia tions ("N owak & Wilmslll999l; 
IWiinan ds. Mille r & van der Kl is 2002). Long term varia- 
tions also have been seen in its optical lightcurve: modula- 
tion over its 9.33 hr orbital period has ranged from ±10% 
and sinusoidal (Thorstensen 1987) to ±30% and complex 
dHakala. Muhh & Dubus. ,1999.) . iHakala. Muhh & DubusI 



1991 " interpret these changes in the optical lightcurve as 



evidence for an accretion disk with a large outer rim, pos- 
sibly due to a warp, being nearly edge-on and partly oc- 
culted by the secondary. The la ck of any X-ray evi- 
dence for binary orbital modulation (|Nowak & Wilmal 19991: 
IWiinands. Miller & van der Klisll2002l) 
bital inclination cannot exceed « 75° 



indicates that the or- 
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Observations with Ginga showed a soft spectrum that could 
be fit with a multi-temperature disk blackbody (diskbb; 
iMitsuda et al.l [1984,) plus a power law tail with photon in- 
dex r « 2-3 jYaqoob etaiT[T99l . The power-law com- 
ponent in those observations comprised up to 25% of the 
1-18 keV flux. Later RXTE observations also were con- 
sistent with a disk blackbody model, but showed no evi- 
dence of either a hard component or any X-ray v ariabil- 
ity above background fluctuatio n s (iNowak & Wilms! [1999). 
IWiinands. Miller & van der KlisI (l2002h conducted a series 
of Target of Opportunity observations designed to catch 
4U 1957+1 1 at the high end of its count rate as determined by 
file All Sky Monitor (ASM) on-board of RXTE. They found 
that at its highest flux, 4U 1957+1 1 exhibited bofli a hard tail 
and mild X-ray variability. 

For all of the above cited X-ray observations, the disk pa- 
rameters had two attributes in common: the best fit inner disk 
temperatures were high (up to 1.7 keV), and the fitted nor- 
malizations were low (~ 10). Nominally, the disk blackbody 
normalization corresponds to (7?km/Dio)^cos6', where 7?km is 
the disk inner radius in km, Dn) is the source distance in units 
of 10 kpc, and 9 is the inclination of the disk. In physical 
models, such large temperatures with such low normalizations 
normally can be achieved only by a combination of large dis- 
tance, high inclination, low black hole mass, high accretion 
rate, and possibly rapid black hole spin. Disk temperatures 
increase as the 1/4 power of fractional Eddington luminosity, 
but decrease as the 1/4 power of black hole mass. Thus, lower 
mass serves to increase the temperature and decrease the nor- 
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malization via decreasing R^m. Large distances and high incli- 
nation also serve to reduce the normalization. (Again, the lack 
of X-ray eclipses limits / < 75°.) Finally, rapid spin yields 
higher efficiency in the accretion disk (and thus higher tem- 
peratures), and a decreased inner disk radius. 

As pointed out by a number of authors, however, the 
diskbb model is not a self-consistent d escription of a 
physical accretion disk (see, for example. jLi et al.l 120051; 
iDaviset al.1 120051; iDavis. Done"&B laes 200^. Several au- 
thors have developed more sophisticated models that incor- 
porate a torque (or lack thereof) on the inner edge of the 
disk, atmospheric effects on the emerging disk radiation field, 
the effects of limb darkening and returning radiation (due 
to gravitational light bending), and the effects of bla ck hole 
spin. Two of these models are the kerrbb m odel of Li et al.l 
(|2005) and the bhspec model of Davis et al. (2005). Using 
these models on soft-state BHC spectra with minimal hard 
tails, various claims have been made as to their ability to 
obser v ationally constrain black ho l e spin (e.g., [Sh afee et al. 
2006; Davis, Done & Blaes| l2006t iMcChntock et al.1 120061: 
Middleton et al. 200C" Given the apparent dominance of a 
soft spectrum in 4U 1957H-11, and the typical weakness of 
any hard tail, 4U 1957H-11 becomes an excellent testbed for 
assessing the ability of these models to constrain physical pa- 
rameters. 

The outline of this paper is as follows. First, we describe 
the analysis procedure for the Chandra, XMM-Newton, and 
RXTE data (©. We then discuss fits to the Chandra data 
with both phenomenological (i.e., diskbb) and more phys- 
ically motivated (i.e., kerrbb) models (©. We then dis- 
cuss the XMM-Newton data, and specifically look at the ab- 
sorption edge regions in both the Chandra gratings spectra 
and t he X MM-Reflection Gratings Spectrometer (RGS) spec- 
tra 03.21 1. We then discu ss sp ectral fits to 48 observations 
from the RXTE archives 04.1l i. and we consider the variabil- 
ity properties of these observations 04. 2b . Finally, we present 
our conclusions and summary (©. 

2. DATA PREPARATION 
2.1. Chandra 

4U 1957H-11 was observed by Chandra on 2004 Sept. 7 
(ObsID 4552) for 67 ksec (two binary orbita l periods). The 
High Energy Transmission Gratings (HETG; ICanizares et al.l 
12005) were inserted. The HETG is comprised of the High 
Energy Gratings (HEG), with coverage from « 0.7-8 keV, and 
the Medium Energy Gratings (MEG), with coverage from w 
0.4-8 ke V. The data readout mode was Timed Exposure-Faint. 
To minimize pileup in the gratings spectra, a 1/2 subarray was 
applied to the CCDs, and the aimpoint of the gratings was 
placed closer to the CCD readout. This configuration reduces 
the readout time to 1 .741 sec, without any loss of the dispersed 
spectrum. 

We used CIAO v3 . 3 and CALDB v3 . 2 . to extract 
the data and create the spectral response files ^ The loca- 
tion of the center of the 0* order image was determined us- 
ing the findzo.sl routine^, which provides « 0.1 pixel 
(« 0.001 A, for MEG) accuracy. The data were reprocessed 
with pixel randomization turned off, but pha randomization 
left on. We applied the standard grade and bad pixel file fil- 
ters, but we did not destreak the data. (We find that for sources 

' We did, however, use a pre-release version of the Order Sorting 
and Integrated Probability (OSIP) file from CALDB v3. 3.0. 
^ |http : / /space .mit ■ edu/ASC/analysis/f indzo/| 



as bright as 4U 1957H-1 1, destreaking the data can remove real 
photon events, while order sorting is akeady very efficient at 
removing streak events.) 

Although our instrumental set up was designed to minimize 
pileup, it is still present in both the MEG and HEG spectra. 
Pileup serves to reduce the peak of the MEG spectra, which 
occurs near 2keV, by approximately 13%, while it reduces 
the peak of the HEG spectra by approximately 3%. In the Ap- 
pendix we describe how we incorporate the effects of pileup 
in our model fits to the Chandra data. 

All analyses, figures, and tables presented in this 
work were produced using the Interac tive Spectral 
Inte rpretation System (ISIS; iHouck & Denicolal 
I2OOOI) . (Note: all 'unfolded spectra' shown in the figures were 
generated using the model-independent definition provided by 
ISIS; see Nowak et al. 2005 for further details.) For our x^- 
minimization fits to Chandra data, we take the statistical vari- 
ance to be the predicted model counts, as opposed to the ob- 
served data count s. (C hoosing the former aids in our Bayesian 
line searchs; see ^3.21 ) Additionally, we evaluated the model 
for each individual gratings arm, but calculated the fit statistic 
using the ISIS combine_datasets function to add the 
data from the gratings arms-^. The Chandra data show very 
uniform count rates and spectral colors over the course of the 
observation, therefore throughout we use data from the entire 
67 ksec observation. 

2.2. XMM-Newton 

4U 1957H-11 was observed by XMM-Newton on 2004 
Oct. 16 (ObsID 206320101) for 45 ksec. XMM-Newton car- 
ries three different i nstruments, the European Photon Imag- 
ing Cameras (EPIC; IStriider et al.llIOOll iTurn er et al.'l2001j) 
the R eflection Grating Spectrometers (RGS; [d en Herder et aT| 
I2001h . and the Optical Monitor (OM; lMason et al.i2001.) . The 
OM was not used in this analysis. 

The EPIC instruments consist of 3 CCD cameras, MOS- 
1, MOS-2, and pn, each of which provides imaging, spectral 
and timing data. The pn and MOS-1 cameras were run in 
timing mode which provides high time resolution event in- 
formation by sacrificing 1 -dimension in positional informa- 
tion. The MOS-2 camera was run in full-frame mode. The 
EPIC instruments provide good spectral resolution (A£'=50- 
200 eV FWHM) over the 0.3-12.0 keV range. There are two 
RGS detectors onboard XMM-Newton which provide high- 
resolution spectra (A/AA=100-500 FWHM) over the 5-38 A 
range. The grating spectra are imaged onto CCD cameras 
similar to the EPIC-MOS cameras which allows for order 
sorting of the high-resolution spectra. 

The XMM-Newton data were reduced using SAS 
version 7.0. Standard filters were applied to all XMM- 
Newton data. The EPIC-pn data were reduced using the pro- 
cedure ep chain. We reduced the EPIC-MOS data using 
emchain. Source and background spectra for the pn and 
MOS-1 data were extracted using filters in the one spatial co- 
ordinate, RAWX. Response files were created using the SAS 
tools rmf gen and arf gen. The MOS-2 data were found to 
be considerably piled up and were therefore not used in this 
analysis. 

^ Owing to the fact that the pileup correction is a convolution model that 
requires knowledge of the individual response functions for each gratings 
arm. the dataset combination performed here is distinct from adding the PHA 
and response files before fitting. (As a result, this combined data analysis is, 
in fact, impossible to reproduce in XSPEC.) 
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Fig. 1 . — Fits to the summed HEG (hollow points) and MEG (solid points) 
data in the 1.8-13 Arange. The neutral column was fixed to be the same for 
all data sets, but pileup connection was treated on an individual detector ba- 
sis, and the HEG diskbb parameters were allowed to differ from the MEG 
di s kbb parameters. 

The RGS data were reduced using rgsproc, which pro- 
duced standard first order source and background spectra and 
response files for both detectors. 

2.3. RXTE 

During the past ten years, RXTE has observed 4U 1957H-1 1 
n umer ous times. Th e first o f these observations was presented 
in iNowak & Wilmsl ( Il999i) . A series of observations, specif- 
ically triggered to observe high flux sta tes, were presented in 
IWijnands. MiUer & van der KlisI (l2002h . One RXTE observa- 
tion was scheduled to occur simultaneously with the Chandra 
observation (PI: Nowak), while another was scheduled to oc- 
cur simultaneously with the XMM-Newton observation (PI: 
Homan). The remaining observations come from different 
monitoring campaigns'*, conducted by various groups. We ob- 
tained all of these data from the archives - 108 ObslDs, 4 of 
which we exclude as they occurred during a brief period when 
RXTE experienced poor attitude control. We combine Ob- 
slDs wherein the observations were performed within a few 
days of one another and color-intensity diagrams indicate lit- 
tle or no evolution of the source spectra, yielding 48 spectra. 

The data were prepared with the tools from the HEA- 
SOFT v6.0 package. We used standard filtering criteria for 
data from the Proportional Counter Array (PCA). Specifically, 
we excluded data from within 30 minutes of South Atlantic 
Anomaly (SAA) passage, from whenever the target elevation 
above the limb of the earth was less than 10°, and from when- 
ever the electron ratio (a measure of the charged particle back- 
ground) was greater than 0.15. Owing to the very modest flux 
of 4U 1957-1- 11, we used the background models appropriate 
to faint data. We also extracted data from the High Energy 
X-ray Timing Experiment (HEXTE), but as 4U 1957H-11 is 
both faint and very soft, none of these data were of sufficient 
quality to use for further analysis. 

We applied 0.5% systematic errors to all PCA channels, 
added in quadrature to the errors calculated from the data 
count rate. For all PCA fits, we grouped the data, starting at 
> 3 keV, with the criteria that the signal-to-noise (after back- 
ground subtraction, but excluding systematic errors) in each 
bin had to be > 4.5. We then only considered data for which 

^ Some of these observations were scheduled simultaneously with ground- 
based optical observations. No optical data are included in this work. 




o : i I I n M I ■ 

0.5 1 2 5 

Energy (keV) 

Fig. 2. — A disk plus Comptonization model (see text) fit to the combined 
MEG data. Dashed lines show the individual model components. 

the lower boundary of the energy bin was > 3 keV, and the 
upper boundary of the energy bin was < 18 keV. This last cri- 
terion yielded upper cutoffs ranging from 12-18 keV. 

3. SPECTRA VIEWED WITH CHANDRA AND XMM-NEWTON 
3.1. Broadband Spectra 

Our main goals in modeling the Chandra data are to de- 
scribe the broad band continuum, accurately fit the edge 
structure due to interstellar and/or local-system absorption, 
and to search for narrow emission and absorption line fea- 
tures. To accurately describe the absorption due to the in- 
terstellar medium, w e use tbnew, an updated versio n of the 
absorption model of I Wilms. AUen & McCravl (l2000h . which 
models the Fe L2 and L3 edges, and includes narrow reso- 
nance line structure in the Ne and O edges, as have been 
observed at high spectral resoluti on with C handra-HETG 
(Juett, Schulz & Chakrabarty 2004; lJuett et al. 2006, abun- 
dances and depletion factors used here have been set to be 
consistent with these previous studies, i.e., they take on the 
default parameter values of the model, with the exception of 
the Fe abundance parameter being set to 0.647, and the Fe and 
O depletion parameters both being set to 1). 

We began by jointly modeling both the HEG and the MEG 
spectra; however, as seen in Fig. [T] we find systematic dif- 
ferences (up to K, 5%, which is less than the statistical noise 
level shown in the figure) between the gratings arms, espe- 
cially at low energies, which cannot be accounted for via 
the pileup modeling. Specifically, the HEG spectra show a 
greater soft excess, such that fits to solely those data in the 
0.7-8 keV range with an absorbed diskbb model yield no 
absorption, contrary to the edge structure detected by MEG 
below 0.7 keV. Given the disagreements with MEG, the lack 
of a fittable column in HEG, the lack of apparent lines in the 
1-8 keV region (see below), and the fact that the MEG allows 
us to model line and edge structure in the 0.5-1 keV region, 
we shall only consider MEG data. 

Systematic disagreements were also found between the si- 
multaneous Chandra and RXTE data. The RXTE data were 
slightly softer in the regions of overlap (k, 4-7 keV), and ow- 
ing to the very large effective area of RXTE, we found that 
the latter instrument dominated any joint fits. The two in- 
struments fit qualitatively similar models with comparable pa- 
rameters; however, no joint model, even when applying sys- 
tematic errors to the PCA data, produced formally acceptable 
fits. The differences between the instruments are likely at- 
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Fig. 3. — Left: The kerrbb model, with a* » 1. fit to the combined MEG data. System distance, inclination, and mass were fixed (see text), but spectral 
hardening factor and black hole spin were left free. Right: same data and model as on the left, except that now the black hole spin is fixed to a* = 0. 



tributable to calibration uncertainties for both, and therefore 
we chose to analyze the Chandra and RXTE data indepen- 
dently. 

For our subsequent broadband spectral fits to the MEG 
data, we grouped the il'*' orders to a combined signal-to- 
noise of > 6 and a minimum of two wavelength channels 
(i.e., 0.01 A) per bin. We then fit the data in the 0.45- 
7.75 keV region. Results are presented in Table [T] Over- 
all, the spectra are well-described by a very mildly absorbed 
(Nh ~ 1-1.2 X 10^' cm~^) diskbb spectrum, with a fairly 
high temperature (kTin = 1.75 keV), but low normalization 
(Adbb = 7.2), with little need for additional features such as 
a narrow or broad Fe Ka line (cf. ^3.2b . Formally, if we add 
a line with fixed energy and width of 6.4 keV and 0.5 keV, 
respectively, the 90% confidence upper limit the line flux is 
7 X 10"^ phcm"^ s"'. 

The low disk model normalization could correspond to a 
combination of low black hole mass, high inclination, and 
large distance. The fitted temperature is on the high side for a 
'soft state' BHC. Fits to soft states of LMC X-1 and LMC X- 
3, for examp le, have consisten tly fo und kTin ^ 1.25 keV 
(lNo\vak et al.ll2001l: IWilms et alj|2001l: iDavis. Done & BlaesI 
|2006|), while our own fits to RXTE and ASCA data of 
4U 1957H-11 have found kTin ^ 1.6keV. Such a high tem- 
perature might be taken as a hallmark of some combination 
of high accretion rate, low black hole mass, and/or high black 
hole spin. Another possibility is the presence of a Comptoniz- 
ing corona, serving to harden any disk spectrum. 

To show, at least qualitatively, how the presence of a corona 
would affect our disk fit parameter s , we a dded the comptt 
Comptonization model ( TitarchukI 1 19941) . to the diskbb 
model, even though the disk fit residuals do not indicate the 
necessity of an extra component. Given the quality of the 
simple disk fit, we chose to reduce the number of free param- 
eters by tying the temperature of the seed photons input to the 
corona to the temperature of the disk inner edge. Furthermore, 
we froze the temperature of the corona itself to 50 keV. The 
additional fit parameters were then the coronal plasma optical 
depth, Tp (with its lower limit set to 0.01), and the normaliza- 
tion of the coronal spectrum. Results of this fit are shown in 
Fig.|2]and are presented in Table [T] 

Although the improvement to the fit is only very mild 
(Ax^ = 1 1), we see that the temperature of the diskbb com- 
ponent drops to 1 .25 ke V, while its normalization constant in- 



creases by more than a factor of 2. The comptt component 
now comprises sa 30% of the implied bolometric flux. Fur- 
thermore, at energies ^ 3 keV it is consistent with a F « 2 
photon index power law; however, this latter fact is tempered 
by the limited bandpass of Chandra at these energies. In a 
"physical interpretation" of these results, we see that a very 
weak, mild corona can qualitatively replace the need for an 
unusually hot, low normalization disk. 

As an alternative for explaining the high disk temperature, 
we can consider spin of the black hole. Although the simple, 
two parameter diskbb model does an extremely good job 
of modeling the MEG spectra, it has a number of shortcom- 
ings f rom a theoretical perspective (see lDavis. Done & BlaesI 
l2006l) . Its temperature profile is too peaked at the inner 
edge, its implied radiative efficiency does not match realis- 
tic expectations, it incorporates no atmospheric physics, and 
it includes no (special or general) relativistic effects. In re- 
sponse to these model shortcomings, several authors have 
developed more sophisticated disk models, specifically the 
k errbb model o f Li et al. (2005) and the bhspec model 
of lDavis. Done & BlaesI (I200Q) . Here we consider the former 
model. 

The kerrbb model increases the number of disk param- 
eters from two to seven, and it allows other effects to be in- 
cluded, such as limb darkening and returning radiation to the 
disk from gravitational light bending (both turned on for fits 
described in this paper). The fit parameters are the mass, spin 
(a*), and distance to the black hole, the mass accretion rate 
through the disk (M), the inclination of the accretion disk to 
the line of sight, a torque applied to the inner edge of the disk 
(here set to zero), and a spectral hardening factor (hd). The lat- 
ter is to absorb uncertainties of the disk atmospheric physics, 
and represents the ratio of the disk's color temperature to its 
effectiv e temperature. Preferred values have been hj ^1.7 
(iLi et a l. 2005; Shafeeetal. 2006; Mc Clintock eTaP 12006). 
while other models (e.g., bhspec) attempt to calculate color 
corrections from first principles (Davi s, Done & Blaes 2006). 

To reduce the fit complexity of the kerrbb model to 
match the two parameters of the diskbb model, one typ- 
ically invokes knowledge obtained from other observations 
(i.e., for system mass, distance, and inclination) and from 
theoretical expectations (i.e., for inner edge torque and hd). 
The former is lacking for 4U 1957-1-11. Here we fix the in- 
clination to 75°, to be consistent with the both the optical 
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TABLE 1 
Broadband Fits to MEG Data 



Nh 
(lO^l 


Adbb 


kTi„ 
(keV) 




M 

(10" gs-i) 


hd 


Pfrac 


Acomptt 

(10--) 


Tp kTo 
(keV) 




1 04^°'^ 

^•"^-0.03 

1 06+<'"2 

^•""'-0.04 

^•"^-0.02 

^•"'-0.05 
1 1O+0.04 
^•^•'-0.04 


7 IS+0.18 
17.87:i;« 


1 74+0.02 
'•'^-0.01 

1 9ir+0.03 

^•^-^-0.03 


n QQQQ+O.OOOO 
"■^'^^'^-0.0002 

"•^-^ -0.003 




n s4»+o.ooi 

u. JH-0_Q 
^•^"-0.03 
'•^°-0.04 


1 9')n+0002 
'•^^"-0.003 

J. 7 

T 41+0.03 
^■^'-0.02 


n 10 0+0.007 

0-1"^'*-0.008 
n 124+0.002 
"•^^^-0.008 
f) 1 ,2+0.002 
-0.007 

-0.008 

126^°°'' 


Q 147+0.001 

"■^^'-0.002 


nni+O.OI 5^ 

"•"^-0.00 -"^ 


1861/1588 
1850/1586 
1862/1587 
1933/1587 
1892/1588 



Note. — The statistic was calculated by using the predicted model counts as the variance. Error bars are 90% confidence values (Ax^ = 2.7) 
for one interesting parameter Parameters in italics were held frozen at that value, diskbb normalization, A^bb. is [(Ri„/1 km)(D/10kpc)]^ cos0. 
For further details of the fits, see the text. 
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Fig. 4. — Flux-corrected XMM-Newton EPIC-pn (squares) and -MOS-1 
spectra (circles), fit with independent absorbed diskbb+diskbb models in 
the 0.7-8 keV range. Dashed lines show the individual absorbed disk com- 
ponents for the fit. 0.05% systematic errors were added to the pn data, while 
1 % systematic en'ors were added to the MOS data. 

and X-ray lightcurve be havior dHakala. Muhli & Dubusll999t 
iNowak & Wilmsl Il999h . Furthermore, we fix the mass to 
3M0 and the distance to lOkpc. We discuss these choices 
further in ^4.11 and ^ but here we note that these latter 
choices allow the closest possible distance to 4U 1957H-11, 
while still maintaining a minimum inferred bolometric lumi- 
nosity of > 3% Ledd- Finally, we set the torque parameter to 0, 
as its value is essentially subsumed by degeneracies with the 
fitted mass accretion rate (Li et al. 2005; Shafee et al. 2006). 

Allowing the spectral hardening factor, hd to remain free, 
we obtain a kerrbb fit to the MEG spectra that is of com- 
parable quality to the simple diskbb fit. We note, however, 
that the fitted value of hd = 1 .22 is sm aller than the commonly 
preferred value of 1.7 (see lDone & Davis 2008), and that the 
black hole spin is found to be near unity (a* = 0.9999). Again, 
the latter is being driven by the fact that the fitted diskbb 
temperature is itself rather high for a soft state BHC system. 
Freezing hd = 1.7, we obtain a slightly worse, but still rea- 
sonable, fit (Ax^ = 71), with a large inferred black hole spin 
(a* = 0.95). If instead we freeze the spin a* = 0, but allow 
the spectral hardening factor to be free, we obtain a in- 
termediate to those of the previous two fits (Ax^ = 30 from 
the best kerrbb fit), but with hd = 3.41. This rather large 
valu e is well-outsi de the normal range of variation (hd =1.4- 
1.8;|Done & Davis 2008). As shown in Fig. [3] however, these 
fits are virtually indistinguishable from one another 

We next consider broadband fits to XMM-Newton data. 



These XMM-Newton observations occured only 6 weeks after 
the Chandra observations, with apparently little source varia- 
tion as measured by the RXTE-ASM in the intervening pe- 
riod. As we further discuss in ^ the XMM-Newton spectra 
showed a flux and spectral hardness comparable to the Chan- 
dra spectra. Thus, we expect the XMM-Newton spectra to be 
qualitatively and quantitatively similar to the Chandra spec- 
tra. We have found that in practice, the XMM-Newton spectra 
are somewhat difficult to fit over their entire spectral range, 
and furthermore, the EPlC-pn and -MOS data do not com- 
pletely agree with one another. For a variety of models that 
we tried, strong, sharply peaked residuals occurred for both 
detectors near 0.5 keV (the O K edge), although these resid- 
uals were most pronounced in the pn data. Above 8 keV, the 
two detectors also show (different from each other, and from 
RXTE) spectral deviations from simple models. Accordingly, 
we only consider data in the 0.7-8 keV range. Additionally, 
we add uniform 0.5% (pn) and 1 % (MOS) systematic errors to 
each spectrum. (These levels were required to obtain reduced 
~ 1-2 with 'simple' spectral models.) 

The only simple models that we found to adequately de- 
scribe the data required an additional spectral component be- 
low 2 ke V, which here we model as a second di s kbb compo- 
nent with kTin « 0.3 keV and Ajbb » 1300 (MOS) and w 2600 
(pn). The required Nh « 3 x 10^' cm~^ for both detectors was 
somewhat larger than that found for MEG. The dominant disk 
component, primarily required to fit the 2-8 keV spectra, had 
comparable parameters to the diskbb fit to the MEG data; 
namely, Adbb ~ 8 and kTi,, « 1 .6 keV. 

Significant residuals clearly remain in the 0.7-2 keV range 
for both the pn and MOS data, and the pn and MOS 
spectra quantitatively disagree with each other Given the 
fact tha t the MEG spectra, as well as prior ASCA spectra 
dNgwak & Wilms 1999) were well described with a simple 
absorbed diskbb model, with lower Nh, we are inclined to 
ascribe a large fraction of the additional required « 0.3 keV 
disk component and the low energy residuals to systematic 
errors in the XMM-Newton calibration. 

With the above caveats in mind, the XMM-Newton data, 
however, confirm that in the 2-8 keV range, the characteris- 
tic disk temperature for 4U 1957H-11 is indeed rather large. 
Furthermore, the XMM-Newton spectra show no evidence for 
any narrow, or moderately broad, Fe Ka line. Similar to the 
Chandra data, adding a 6.4 keV line with 0.5 keV width (both 
fixed) yields a 90% confidence upper limit to the line flux of 
6 X 10"^phcm"^ s"'. 

3.2. Edge and line fits 
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TABLE 2 
Absorption Lines Fit to MEG and 
RGS Data 



Line 


A 

(A) 


Eq. Width 
(mA) 


NelX" 
Nelll 
Nell 


1 1 44Q+0.005 
'-'■^^^-0.005 

14.508 


'■-'-3.1 


OVIII K/3'' 
OVIII Ka 
Ovil 
? 

On 


18.629 
18.967 
21.602 
21.80^;«^ 
23.375 


_Q Q+9.9 
^■^-9.0 

-18!1« 
-71+36 
-15+15 

'-^-18 


OVIII KP" 

UVill rs-Cx 

Ovii 
? 

On 


18.629 

If? 0<<7 
J 0. yO/ 

21.602 
23.375 


0+0 

70+15 
2S_i4 

-36+2' 

JO 24 
-67+^ 


Ovill K/3'' 
OVIII Ka 
Ovil 
? 

On 


18.62t«:^ 
18.96^;i 
21 59+«"3 

91 Sl+0.03 
23.^75 


-8.5!|5 
-14!'" 
-19+" 
-23+" 
-29+" 


Note. — Italicized parameters were 
lield fixed. Errors are 90% confidence 
level. 

^ MEG data fit between 13A-15A, 
grouped to S/N> 5 per bin.'' MEG data 
fit between 18A-25A, grouped to 16 
channels per bin.*-' MEG data fit be- 
tween 1 8 A-25 A, grouped to 32 chan- 



nels per bin.'' RGS data fit between 18 A- 
25 A (RGSl) and 18 A- 19.9 A (RGS2), 
grouped to S/N> 5 and > 4 channels per 
bin. 

We have employed a variety of techniques to search for 
Unes in the 4U 1957+1 1 data, ranging from direct fits of nar- 
row components at known line locations to a "blind search" 
using a Bayesian Blocks technique. The latter involves fitting 
a continuum model to the binned data, using the predicted 
model counts as the statistical variance in the fits, and then un- 
binning the data and comparing the likelihood of the observed 
counts to the predicted counts. (For a successful application 
of this technique to Chandra-HETG data of the low lurn inos- 
ity active galactic nuclei, M81*, see lYoung et al.ll2007l ) No 
candidate features that appear to be local to the 4U 1957-1-1 1 
system were found. This is in contrast to, for example, Chan- 
dra observations of GROJ 1655-40, which at a similar source 
luminosity in a spectrally soft, disk-dominated state showed 
very str ong absorption fea tures associated with the disk atmo- 
sphere (iMiller et al.ll2006h . Comparable equivalent width fea- 
tures would have been very easily detected in our observation 
of 4U 1957H-11. We note, however, that GROJ1655-40 itself 
has not always shown disk atmosphere absorpt ion features. 
Two w eeks prior to the observation described by lMiller et alJ 
(I2OO6I) . despite the source being only slightly fainter while 
also being in a soft, disk-dominated state, Chandra-HETG 
observation s revealed no narrow emission or absorption lines 
(iMiller et al.,,2006. .2008.) . The 'duty cycle' with which such 
lines are prominent in disk-dominated BHC spectra compa- 




13.5 14 14.5 

Wavelength (A) 



Fig. 5. — Flux corrected spectra (top panel) and count rate per angstrom 
(middle panel) for the combined MEG ±1'" order spectra, with an absorbed 
diskbb model fit in the 13-15 A region (see text). The top and middle 
panels show the data fit with the Ne II, Ne III, and Ne IX lines, while the 
residuals show a fit without the lines. 

rable to the specti-a of both GROJ 1655-40 and 4U 1957H-1 1, 
remains an open observational question. 

On the other hand, the above Bayesian Blocks procedure 
yielded a few candidate features in the Chandra data that are 
likely attributable to absorption by the Interstellar Medium 
(ISM). The strongest such feature is the absorption Une at 
13.44 A. This feature is consistent with Ne IX absorption from 
the hot phase of the interstellar medium (Juettet al. 2006). 
Accordingly, we also simultaneously fit expected features 
from Ne II and Ne III absorption at 14.61 A and 14.51 A, re- 
spectively. 

We again choose an absorbed disk model, albeit with the 
pileup fraction fixed to the values from our broad-band mod- 
els, and only consider the 13-15 A region. We group the spec- 
trum to a minimum signal-to-noise of 5 in each bin, and we 
let the equivalent widths and wavelengths of the Ne IX and 
Ne II lines be free parameters while fixing the wavelength of 
the Ne III absorption to 14.508 A. Results from this fit are 
presented in Fig. |5] and Table |2] The fitted wavelengths of 
the Ne IX a nd Ne II lines agree well with the MEG studies of 
lJuett et all (f2006), although we find the Ne IX line wavelength 
to be redshifted by 0.005 A (i.e., 1 1 1 km s"'). Our fits, how- 
ever, are within 0.002 A of the theoretical value of 14.447 A. 

Formally, the fit improves with addition of any of the three 
Ne lines, although the Ne III region residual is far from the 
strongest in the overall spectrum. The Ne IX and Ne II ab- 
sorption lines are more clearly detected. Their equivalent 
widths are near the mid-level to higher end^ of the sample dis- 
cussed by Juett et al. (2006). As we discuss in greater detail 
elsewhere (Yao et al.ll2008l) . the 4U 1957H-11 Ne IX equiva- 
lent width is consistent with the source's probable location 
outside the disk of the Galaxy. Specifically, presuming that 
4U 1957-1-11 is greater than 5kpc distant, and choosing a 
model where the hot phase of the ISM lies in layers predom- 
inantly outside the galactic plane, the equivalent width of the 
Ne IX line in 4U 1957-1-11 is consistent with the equivalent 
width of the same feature observed in the Active Galactic Nu- 
clei Mkn 42 1 , appropriately scaled for that sightline through 

' The Ne IX line equivalent width measured in 4U 1957+11 is a fac- 
tor of two l ess than that observed in GX 339—4; however, as discussed by 
IMiller et aP fTOOA) and juett et aL i2006) the latter source's Ne IX column is 
likely dominated by a warm absorber intrinsic to the binary. 
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Fig. 6.— Flux corrected RGS spectra in the 18-25 A (RGSl) and 18- 
19.9 A (RGS2) range, fit with a model consisting of interstellar absoiption, 
a phenomenological continuum model (a broken power law plus multicolor 
blackbody), and various absorption lines (see Table |2). The residuals show 
the fit without the lines. 

our Galaxy. Therefore it is likely that 4U 1957+1 1 is sampling 
a large fraction of the hot phase interstellar gas that lies within 
and directly outside the galactic plane ( Yao & Wang 2007). 

Consistent with this picture, both the RGS and MEG spec- 
tra show evidence of longer wavelength features that we also 
associat e with absorp tion by the interstellar medium (Table|2l 
see also lYao et al.ll200 8). For the RGS data, we grouped both 
the RGSl (fit between 18-25 A) and RGS2 (fit between 18- 
19.9 A) data to have a signal-to-noise of 5 and a minimum of 
four channels per bin. The MEG data were grouped to ei- 
ther 16 channels per bin, or 32 channels per bin. The MEG 
data were fit with an absorbed disk blackbody, while the RGS 
data required additional continuum structure (which we mod- 
eled here with a broken power-law). For all these model fits 
we added absorption lines expected from the warm and hot 
phases of the interstellar medium. The O II line from the 
warm phas e of the ISM (not included in the tbnew model; 
Huett. S chulz & Chakrabartvll2004) is clearly detected by the 
RGS data, and consistent limits are found with the MEG data. 

The hot phase of the ISM is detected with both the RGS and 
MEG via O VIII Ka absorption (18.96 A) measurements. The 
presence of O VIII Kf3 (18.63 A) is not strictly required, but 
the limits on it s equivalent wid th are consistent with the mea- 
sured Ka line (I Yao et al.ll2008l) . In addition to O VIII, we also 
find evidence for O VII Ka absorption. Interestingly, an even 
stronger feature is detected at 21.8 A. We have no good iden- 
tification for this feature (it is 3000 km s"' redshifted from 
O VII). Its equivalent width might change between the MEG 
and RGS observations, and therefore it could be a candidate 
for a line intrinsic to the source. 

It is interesting to note that for the 18-25 A fits discussed 
above, we find that the implied neutral column is somewhat 
larger than that found with the broadband fits. Specifically, we 
find Nh = 1.9± 0.2, 1 .8 ± 0.5, and 2.0 ± 0.2 x 10-' cm'^ for 
the RGS and MEG (16 channels per bin, 32 channels per bin) 
data. These values are primarily driven by the O edge region 
of the fits, although we have found for our broadband fits that 
merely changing the O abundance is insufficient to bring their 
fitted neutral columns into agreement with the 18-25 A fits. 
It is possible that an additional very low energy component 
added to the broad band fits might allow for a larger column 
that is consistent with these narrow band fits. 
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Fig. 7. — Flux corrected MEG spectra in the 18-25 A range, fit with a 
model consisting of interstellar absorption, a multicolor blackbody, and var- 
ious absorption lines (see Table |2). The residuals show the fit without the 
lines. The top panel shows the spectra grouped to 16 channels per bin, while 
the bottom panel shows the spectra grouped to 32 channels per bin. 

4. RXTE OBSERVATIONS 

In the previous sections we have shown that the 
4U 1957+11 spectrum is consistent with a disk having a 
high peak temperature. One could reduce the need for this 
large peak temperature by including a hard component, e.g., 
a Comptonization spectrum; however, this added component 
only becomes dominant at photon energies above the band- 
passes of Chandra and XMM-Newton. Thus, to distinguish 
between the possibilities of high black hole spin or an addi- 
tional hard spectral component leading to the high fitted tem- 
peratures, we turn toward the RXTE data, which can provide 
spectral information up to « 20keV for4U 1957+11. 

4.1. Spectra 

In Fig. [8] we present the RXTE-ASM lightcurve, with 
6 day bins, and also show the scaled 3-18 keV flux from the 
pointed RXTE observations. This figure highlights the times 
of the Chandra and XMM-Newton observations (which were 
simultaneous with RXTE ObslDs 90123-01-03 and 90063- 
01-01, respectively). Generally speaking, 4U 1957+11 ex- 
hibits an X-ray lightcurve with a factor of 3-A variabil- 
ity on a few hundred day time scale. iNowak & Wilmll 
(Il999h hypothesized that this long term variability was quasi- 
periodic; however, our longer lightcurve does not show 
any clear super-orbital periodicities. Note that the Chandra 
and XMM-Newton observations occur during the same lo- 
cal peak in the ASM lightcurve and have very similar 3- 
ISkeV fluxes as measured by the PC A. We shall present 
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TABLE 3 

(diskbb+gaussian) +powerlaw, OR +comptt FITS TO SELECTED RXTE DATA 



ID 


MJD 


Adbb 


T 


Api 


r 


Acomptt 

(10-^) 




■^line 










(keV) 


(10-2) 






(10-») 


(keV) 



40044-01-02-01 


51271 


5 05+0-'*' 
-0.47 


1 on 1+0.022 


28.54!5.45 


9 1 ^+0.07 
^•"-0.07 






3 9+'-^ 


0+0.6 
-0.0 


48.9/28 






24.98+0.02 


1 226+0001 






249 8+0.8 


A 99+0.00 
"•^^-0.00 


2 4+1.2 


o.o^:« 


25.1/28 


70054-01-11-00 


52656 




1 oio+O.OOl 


00+0 45 
-0.00 








3 n+0.5 




16.7/18 








I 1 +0.001 
^•-"^ -0.001 






0+'.' 
"• -0.0 


53+^."" 
-0.52 


4 2+0.5 


0.8:0.2 


16.4/18 


90123-01-03 


53255 


11.49!0.24 


1 ct:4+0.008 


71+4.31 
' ^-0.70 


2.38:0.62 






4 6+19 


5+0.3 
"•■'-0.4 


25.7/26 






12.10!0:03 


1 cao+0.026 






5 6+0.4 


01+4.99 
"•"^-0.00 


4 g+I.8 


5+0.3 


25.3/26 


90063-01-01 


53294 


11.48!0.29 


1 riro+fl.OOV 


1 61+0.02 


2.60:0.40 








6+0.2 

"•"-0.3 


19.0/26 






12.24!0;42 


1 590+0.029 






6.4:0.8 


01+4.99 
"•"^-0.00 




6+0.2 

"•"-0.3 


18.9/26 



Note. — diskbb normalization, Ajbh. is [(Rin/l km)(D/10kpc)]2 cos9. Power law normalization, Api, is photons/cm^/s/keV at 
1 keV. Line normalization, Aii„e, is total photons/cm- /s in the line. 0.5% systematic errors were added in quadrature to the data. Error 



bars are 90% confidence level 




Fig. 8.— All Sky Monitor lightcurve for 4U 1957+11 with scaled 3-18 keV 
PCA flux from pointed observations overlaid. RXTE observations simulta- 
neous with Chandra and XMM-Newton are circled. Circles here and in sub- 
sequent figures represent observations where 4U 1957+11 was transiting to 
'steep power-law state' behavior. 

detailed results for these two observations, as well as for 
the faintest RXTE observation (ObsID 70054-01-11-00) and 
the brightest RXTE observation (ObsID 40044-01-02-01). 
This latter observation has been pre sented previously by 
IWiinands. Miller & van der Kli s' ('2002'). 

As discussed bv iNow ak & W ilms! d 19991) and 
IWiinands. Miller & van der Khsl (i2002), the RXTE data 
of 4U 1957-1-11 can be well-described by a model that 
consists of a combination of a multi-temperature disk 
spectrum, a power-law, and a gaussian line. For most of 
the RXTE observations the line peaks at an amplitude of 
^ 2-3% of the local continuum level, and it is statistically 
required by the data. No such line was present in either 
the Chandra or the XMM-Newton observation (line flux 
^ 7 X 10"^ phcm"^ s"'). There are several possibilities 
for this discrepancy, and perhaps some combination of all 
are at play. With its « 1 deg^ effective field of view, the 
PCA will detect diffuse galactic emission which shows 
a prominent 6.7 keV line. Based upon scalings between 
the PCA fine flux and flie DIRBE 4 .9 Mm infrared flux 
jRevnivtsev. Molkov & Sazonovl 120061) . one expects an 
RXTE-detected line amplitude of « 4 x 10~''phcm"^s"' 
(i.e., approximately 10% of the observed line) at the Galactic 
coordinates of 4U 1957H-11 (I = 51.307°, b = -9.330°). 



Remaining systematic uncertainties in the PCA response 
matrix may be as large as 1 % of the continuum level in the 
Fe line region, which could account for a line amplitude of 
^ 1.5 X lO'^^phcm"^ s"' for the observations simultaneous 
with Chandra and XMM-Newton. The remaining 50% of 
the line flux for these latter two observations, which is still 
three times larger than the Chandra and XMM-Newton upper 
limits, may be real and could be indicative of the greater 
sensitivity of PCA to weak and broad line features. Such 
a broad line would have been detectable by Chandra and 
XMM-Newton only if their own systematic uncertainties 
were at the 1 % level, and we have already noted larger broad 
band discrepancies between the HEG and MEG. 

For the disk plus power-law fits, we constrained the power- 
law photon index, F, to lie between 1-3, but let the power-law 
normalization freely vary. These fits were fairly successful 
overall (reduced ranged from 0.5-2, and averaged 1). Pa- 
rameter trends with observed flux are shown in Fig. |9] and 
show several features of interest. First, the vast majority of 
the observations are consistent with having a nearly constant 
disk normalization, which one could take as a proxy for disk 
radius. At the low end of 3-18keV flux, the disk normaliza- 
tion/radius turns upward. Phenomenologically, such observed 
increases in disk radius have been associated with low lumi- 
nosity transition s to the hard state ( e.g., see the similar behav- 
ior of LMC X-3: Wilms et al.l2001l) . wh ich we expect to occur 
near 3% of the Eddington luminosity (iMacc arone 2003). It 
is therefore possible that our lowest luminosity observations 
are near su ch a transition at 3% LEdd- (On the other hand, 
ISaito et aH 12006 note a similar increase in fitted disk radius 
for GRO J 1655-40, without any associated transition, or near 
transition, to a low/hard state.) 

Second, we see that as the 3-18 keV flux increases, the 
fraction of the flux contributed by the power-law slightly in- 
creases. Half a dozen points, however, show a much more sig- 
nificant contribution by the power-law (Fig. |9ll. At the same 
time, those observations show a dramatic decrease in the disk 
normalization. Each of those observations is associated with 
peaks in the ASM lightcurve, and several are associated with 
high variabi lity, including one instance of "rapid state transi- 
tions" (see ^4.2ll. We ide ntify those observations with what 
iRemillard & McClintockl (2006) ref er to as the "steep pow er- 
law" state (or, the "very high state"; Mivamo to et al.lll993h . 

Even though the disk normalization drops and the power- 
law flux increases for the disk plus power-law fits, the fit- 
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TABLE 4 

kerrbb+comptt+gaussian FITS TO SELECTED RXTE DATA 



ID 


a 


M 

(10"gs-l) 






Dbh 
(kpc) 


A-comptt 

(10-^) 




(10-*) 


""line 

(keV) 


xVdof 


40044-01-02-01 


0.99026!<i['S 


0.749!}i-"' 

-(). 1 1 J 


n 029+0.045 
"■^^^-0.043 


10 


263.7+" 

-8.5 


0.19!«o? 


5.4+4-6 

—3,0 


—0.5 


18.0/27 




0.99409+['!]!!5«i 


1.541+^?'j 

-U. /61 




1.7 


-2.8 


264.0+'?„5 

-10.4 


-0.0 1 


^■^-3. 1 


0.5+0-5 

-0.5 


18.0/27 




0.59987+!'08325 

-0.08022 


3.556+0-471 

-0.452 




1.7 


10 


254.2+"- 

-16.5 


0.20!*!'' 

-0.02 


4.3!|9 


0.3+0-4 

-0.3 


19.0/28 


70054-01-11-00 


0.99896!<'!'S 


0.349!<'-001 

— U.UUl 


1 


1 09+0.002 
^"^-0.003 


10 


— u.u 


0.01+4-95 
— u.uu 


— U.o 


0.8+0-2 

— U.^ 


16.7/17 




99986+'"«"«» 


1 151+0.002 

1.1D1_0 004 




1.7 


19 5+"-' 
^^■■^-0.0 


4 6+0 2 
^■"-0.4 


001+0.01 

"■"^-0.00 


7 1+0.4 
'■^-0.7 


1 0+00 
'•"-0.1 


25.7/17 






1 414+0.035 

^■^^^-0.019 




1.7 


10 




001+4.99 
"■"^-0.00 


4 5+1-7 
^■^-1.0 


7+0-2 
"•'-0.2 


46.5/18 


90123-01-03 


90069+0 00021 


Q 576+0.467 
U.J ;o 017 


1 


no7+0.017 
'-0.014 


10 


5.6^:« 


02+2 50 
"•"-'^-0.01 


4 7+' -9 
^■'-1.5 


5+0-3 
"•^-0.4 


25.0/25 




n 9Q974+0.00012 
-0.00014 


2 676+0- 154 
z. 0/0^016 




1.7 


2 1 7+0.4 
^^■'-0.1 


6 0+<' 7 

D.u_3 


01+4-99 
"■"^-0.00 


^■6U| 


5+0-2 


24.9/25 






y 419+0.033 




1.7 


10 


0+2 1 


37+4.63 
"■^'-0.36 


9.0!1:5 


7+0.1 
"•'-0.1 


48.0/26 


90063-01-01 




655+''-"'*2 

V.VJJ^ O Q21 


1 


1 , 1+0.021 
^ ^-"-0.011 


10 


ft 9+0.8 


01+" 68 

"•"^-0.00 


5 l+'-6 


6+0-2 

"•"-0.3 


18.8/25 




n 9979<;+0.00050 
u.y?;7J_g 1,0268 


1 7nQ+0.474 
i./U8_0 472 




1.7 


23.3!?;2 




03+4-97 
"■'-"-0.02 


5 l+i-'* 


6+0" 

"•"-0.3 


18.8/25 




n 0/1 SI £+0.00203 
U.S451b_(, 00401 


2415+0.035 
^■^'-'-0.026 




1.7 


10 


0+2-3 


2 50+2-^" 


8.7+[;6 


7+0.1 

"•'-0.1 


39.0/26 



Note. — One fit used a fixed mass of 3 Mq , a fixed distance of 10 kpc, and a variable spectral hardening factor, another set used a 
fixed mass of 16 Mq, a fixed spectral hardening factor of 1.7, and a variable distance, and the third set a fixed mass of 3 Mq, a fixed 
distance of 10 kpc, and a fixed spectral hardening factor of 1.7. Line normalization, Aii„j., is total photons/cm2/s in the line. 0.5% 
systematic errors were added in quadrature to the data. Error bars are 90% confidence level. 
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Fig. 9. — Left: The fraction of disk flux to total flux, in the 3- 1 8 keV energy band, for diskbb+powerlaw fits to 4U 1957-1-11. Middle: Square root of disk 
normalization vs. PCA flux for diskbb+powerlaw fits to RXTE observations of 4U 1957-1-1 1. Right: Disk temperature vs. PCA for diskbb+comptt fits 
to4U 1957-1-11. 

3 M0. The latter choices are motivated by our faintest RXTE 
observation for which we measure a 3-18 keV absorbed flux 
of 3.1 X 10"'" ergs sec"' cm"^. Taking the (diskbb model at 
face value, with an inclination of 75°, this result translates to 
an unabsorbed bolometric luminosity of 3% LEdd for this pre- 
sumed mass and distance. The black hole is unlikely to be less 
massive than 3 Mq, and low black hole masses lead to higher 
disk temperatures. Thus, we see these parameter choices as 
minimizing the need for high spin parameters in the kerrbb 
fits, with lOkpc then being a lower limit to the source dis- 
tance. 

The remaining kerrbb fit parameters are the accretion 
rate, M, the black hole spin, a*, and the disk atmosphere 
color-correction factor, hd. The kerrbb model does not 
fit a disk temperature per se; therefore, we set the input 
seed photon temperature of the comptt model equal to 
the diskbb pe ak temperature from our previous fits (see 
iDavis. Done & B laes 2006). As before, we also freeze the 
coronal temperature to 50 ke V, and constrain the optical depth 
to lie between 0.01-5. Selected fit results are presented in Ta- 
ble Hand Fig. [TO] 

Even with choosing a very low black hole mass, and thus 
naturally favoring higher disk temperatures, the kerrbb 
models uniformly prefer to fit a* « 1. On the other hand. 



ted disk temperature does not show any significant deviations 
from a trend of 3-18keV flux cx (kT)^ (not shown in Fig.|9j. 
For a disk spectrum with constant disk radius, we expect the 
bolometric flux to scale cx (kT)*, and, for these temperatures, 
to scale approximately as (kT)^ in the 3-18 keV bandpass. 
This same relation between disk temperature and 3-18 keV 
flux holds if we instead fit the spectra with a disk plus Comp- 
tonization spectrum, as shown in Fig.|9] 

For these latter fits, we model the spectrum in the same 
manner as in |3] we use a diskbb + comptt model and 
tie the seed photon temperature to the disk peak tempera- 
ture, freeze the corona temperature to 50keV, but let the coro- 
nal spectrum normalization and optical depth (constrained by 
0.01 < Tp < 5) be fit parameters. Results for several of these 
fits are presented in Table[3] The majority of the observations 
again show a nearly constant disk radius, and a 3-18 keV flux 
that scales approximately oc (kT)^. The Comptonization mod- 
els, however, now show the half dozen 'steep power-law state' 
observations as being associated with dramatic drops in the 
disk/seed photon temperatures (Fig.|9]l. 

We now consider fits where the unphysical diskbb com- 
ponent is replaced with the kerrbb model. As for the fits 
to the Chandra data, we fix the disk inclination to 75°, and 
set the source distance to 10 kpc and the black hole mass to 
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Fig. 10. — Count rate spectra for the RXTE observations simultaneous with Chandra (top left) and XMM-Newton (top right), and for 
(bottom left) and hardest/brightest (bottom right) of the RXTE observations, all fit with the kerrbb+comptt+gaussian model (3M0, 
varying color correction fraction). Dashed lines show the individual model components folded through the detector response. 

the disk atmosphere spectral hardening remains low, with 
hd w 1.1. There is a slight trend for hd to decrease with in- 
creasing flux, as shown in Fig. [TT| For some of the observa- 
tions where the disk contribution to the total flux dramatically 
decreased, hd drops slightly lower still. This is not surpris- 
ing as the inclusion of a coronal component can qualitatively 
replace the need for hardening from the disk atmosphere. 

The kerrbb model has a number of parameter degenera- 
cies which we can exploit to search for fits with the more usual 
color-correction factor of hd = 1.7. Specifically, the apparent 
temperature increases cx hd, and cx (M/M^)'/'*. If we increase 
hd, then in order to retain roughly the same spectral shape 
we must decrease (M/M^)'/'*. Solely decreasing M, however, 
makes our faintest observation less than 3%LEdd- Addition- 
ally, the source distance would need to be reduced oc M'/^ in 
order to retain the same observed flux. We would have ex- 
pected to detect a transition to a hard spectral state if M were 
lower than we have assumed. In order to keep our faintest 
observation at « 3%LEdd, we must instead scale M cx hd^, 
M oc hd**, and the source distance as cx hd^. A consistent pic- 



the softest/faintest 
lOkpc, and freely 



6x10 8x10 
3- 18 keV Flux (eri 



ture may be obtained with hd = 1.7, if M w 16 M0 and the 
distance to the source is « 22 kpc. 

We searched for a set of kerrbb+comptt fits where we 
froze the disk color-correction factor at hd = 1.7, the black 
hole mass at 16 M0, but let the disk accretion rate, M, spin 
parameter, a*, and source distance all be free parameters. 



Fig. 1 1 . — Spectral hardening factor vs. PCA flux for kerrbb+comptt 
fits to 4U 1957+11. We assume a black hole mass of 3 Mq, a distance of 
lOkpc, and an inclination of 75°. 

Clearly, the source distance cannot truly be physically chang- 
ing in a perceptible way; however, the degree to which we 
obtain a uniform set of fitted distances might provide a self- 
consistency check on the use of the kerrbb model. Selected 
results from these fits are presented in Table [H and the fit- 
ted distances vs. 3-18 keV flux are presented in Fig. [12] In 
general, these fits work every bit as well as the 3 M0, 10 kpc 
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Fig. 12. — Fitted black hole distance for kerrbb+comptt fits to 
4U 1957+11, assuming a spectral hardening factor of hj = 1.7, a black hole 
mass of 16 Mq, and an inclination of 75°. 



fits. The lower flux observations cluster about a fitted dis- 
tance of w 22kpc. There is a trend, however, for the fitted 
distances to decrease with higher fluxes, and for several of the 
'steep power-law state' observations to fit distances as low as 
15kpc. 

We have considered another class of fit degeneracy: disk at- 
mosphere hardening factor and spin. Here we again freeze the 
black hole mass to 3 M© and the distance to lOkpc, but freeze 
the hardening factor to hd = 1 .7. As more of the peak in the ap- 
parent temperature is being attributed to the color-correction 
factor, the disk spin must be decreased, thereby decreasing 
the accretion efficiency and increasing the radius of the emit- 
ting area near the peak temperatures. To maintain the same 
observed flux, the accretion rate must be increased above that 
for the models with lower hd. Formally, these are significantly 
worse fits than our previous ones, with a mean Ax^ = 22. I.e., 
an increased accretion rate and spectral hardening factor are 
not acting as simple proxies for high spin in these models. (In 
terms of fractional residuals, however, the fits are somewhat 
reasonable.) 

As shown in Fig. [13] these fits yield a relatively uniform 
spin of a* « 0.84, with a slight trend for fitted spin to increase 
with increasing flux. Notable exceptions to this behavior are 
seen, however, as some of the 'steep power-law state' obser- 
vations fit much reduced black hole spins. Again, this is a con- 
tribution from a Comptonization component, which is strong 
for these observations, replacing some of the need for spectral 
hardening due to rapid spin. 

Note that we did search for a set of disk atmosphere model 
fits where we kept the source mass low (3Mq) and moved 
the distance further out (17kpc) while fitting higher accre- 
tion rates, to intrinsically increase the disk temperatures and 
raise the implied fractional Eddington luminosities by a fac- 
tor of « 3. These models failed completely, with best fit 
values increasing by factors of almost 3 over high-spin mod- 
els. We again find that high-spin is not simply degenerate 
with increased accretion rate in the models. Essentially this 
is because high-spin adds an inner disk region with high tem- 
perature and high flux (approximately tripling the flux when 
going from non-spinning to Kerr solutions), which leads to a 
strong, broad component in the high energy spectrum. High 
accretion, low spin models, fail to reproduce such broad, high 
energy PCA spectra. 
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Fig. 13. — Fitted black hole spin for kerrbb+comptt fits to 
4U 1957+11, assuming a spectral hardening factor of hj = 1.7, a black hole 
mass of 3 M0 , a distance of 1 kpc, and an inclination of 75 ° . 

4.2. Variability 

The study of the rapid variability properties of 4U 1957+1 1 
was performed using high time resolution data from the PCA. 
We created power spectra from the entire PCA energy band 
(?a2-60 keV), covering the 2"^-2048 Hz frequency range. 
Observations with high disk fractions (>80%, as defined by 
Fig. |9]l typically showed weak or no significant variability. 
Combining these observations resulted in a power spectrum 
that could be fitted reasonably well by a single power-law 
component (x^/dof=101. 0/79). This fit improved significantly 
(X^/dof=73.1/74) by adding two Lorentzians, one (with Q- 
value fixed at Hz) for a weak band-limited noise compo- 
nent around 2 Hz and the other for a QPO at 25 Hz. The re- 
sulting parameters from the fit are listed in Table |5] Note, 
that although the improvement in the fit is statistically signif- 
icant, the two added components are only marginally signifi- 
cant themselves (« 3cr). 

Observations with lower disk fractions showed an increase 
in variability, with maximum of 12.1 ±0.6% root mean square 
(rms) variability (0.01-100 Hz) in the observation with the 
lowest disk fraction (ObslD 50128-01-09-00; disk fraction of 
w30%; see Fig. [I4] |. Combining the power spectra of all ob- 
servations with disk fractions lower than 80% and fitting with 
the same model as above, we find that all variability compo- 
nents increased in strength (see Table |5]). 

The power spectral properties of the high- and low-disk 
fraction observations are consistent with the soft state and the 
soft end of the transition between the hard and soft states, re- 
spectively. Observation of such variability strengthens the ar- 
gument that 4U 1957+1 1 is indeed a black hole, and not a neu- 
tron star. The frequencies of the QPOs are similar to that seen 
in a few other black hole systems when t hey were in or close 
to their soft state, e.g XTE J1550-564 (iHoman et al.ll200l 
and GRO J 1655-40 (ISobczaketal.1 120001). The comparable 
frequencies observed in these latter two sources have ranged 
from 15-22 Hz, and futhermore the spectra of these sources 
have been fit with disk models that imply black hole spins 
in the range of a * « 0.1-0.8 (iDavis. Done & BlaesI 120061: 
IShafee et al.ll2006l) . Thus, if such variability features are re- 
lated to black hole spin, a consistent picture arises between 
their spectral and variability results. 

In two of the low-disk fraction observations (ObslDs 
40044-01-02-01 and 70054-01-04-00) we observed fast 
changes in the count rate, which resemble the dips/flip-flops 
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Fig. 14. — An approximately 25 Hz quasi-periodic oscillation seen in the 
4.5-27.5 keV 4U 1957-1-11 lightcurve from ObsID 50128-01-09-00, the ob- 
servation with the smallest disk contribution to the observed 3-18 keV flux. 
The top panel shows the power spectrum with the Poisson noise level unsub- 
tracted and normalized to 2. The bottom panel shows Fourier frequency 
times the power spectrum, now with the Poisson noise level subtracted and 
the PSD normalized to (root mean square)^/Hz. 



TABLE 5 
Power spectra fit parameters 



Parameter" High Disk Frac. Low Disk Frac. 



PL rms (%) 


1.05±0.09 


2 6+^' 


PL index 


1.57±0.15 


0.92±0.12 


BLN rms (%) 


1.2±0.2 


4 !+'■' 
^■'-0.7 


BLN u„,a, (Hz) 


1 9+' 5 


13!« 


QPO rms (%) 


2.0±0.3 


3.4±0.9 


QPO Q-value 


2 10+' •^'^ 
^■'"-0.76 


2 5+3.6 


QPO v,„c„ (Hz) 


25.0±2.2 


29.5±2.1 


total rms'' 


2.7±0.4 


6.1±0.4 


X^/dof 


73.1/74 


65.9/80 



^ PL = power-law; BLN = band-limited noise; QPO = quasi- 
periodic oscillation'' O.OI-IOOHz 



observed in other black hole X-ra y binaries during transi tion 
from the hard to the soft state ('e.g.. lMiyamoto et alJll991h . In 
4U 1957-1-11, these changes in the count rate were accompa- 
nied by moderate changes in the spectral variability proper- 
ties, but count rates were too low to classify the power spectra 
in the high and low count rate phases independently. 

5. DISCUSSION 

We have presented observations of the black hole candidate 
4U 1957H-11 performed with the Chandra, XMM-Newton, 
and RXTE X-ray observatories. All of these observations 
point toward a relatively simple and soft spectrum, indica- 
tive of a classic BHC disk-dominated soft state. The Chan- 
dra and XMM-Newton spectra, especially at high resolution, 
indicate a remarkably unadulterated disk spectrum. That is, 
the absorption of the spectrum is very low at only Nh = 1- 
2 X 10^^ cm"^, and with the possible exception of an uniden- 
tified line at 21.8 A, there is very little evidence of spectral 
complexity intrinsic to the source. From that vantage point, 
4U 1957-1-1 1 may be the cleanest disk spectrum with which to 
study modern disk atmosphere models. 

Perhaps the one indication of additional, unmodeled spec- 



tral complexity is the fact that the broadband fits and the more 
localized, high resolution edge and line fits yield (low) neu- 
tral columns that differ from each other by a factor of two. 
This, along with the 21.8 A absorption feature might indicate 
the presence of other weak, unmodeled spectral components 
at very soft X-ray energies. Even if this were the case, this 
is to be compared to other BHC to which such atmosphere 
mode ls have been applied, e.g., GRO J1655-40 ( Shafe e et alJ 
120061) . where there is both a larger neutral column (Nh = 7 x 
10^' cm"^), a complex Fe K g region (mod eled with both emis- 
sion and smeared edges by IShafee et al.L2006) and high res- 
olution observations of intermittently present, co mplex X-ray 
spectral features that are local to the system (i.e.. lMiller et alJ 
|2()06, 2008). In comparison, the 4U 1957-1-11 spectrum, even 
with its uncertainties, is much simpler 

The observations presented here strengthen the case that 
4U 1957-1-11 is a black hole system. The spectrum closely 
matches a classic, soft state disk spectrum (well-modeled 
by diskbb), with a contribution from a harder, non-disk 
component that generally increases as the flux increases. 
Occasionally (a half dozen observations), the contribution 
from this hard component increases, the X-ray variability in- 
creases, a possible QPO is observed, and we even observe 
lightcurve 'dipping' and 'flip -flop' behavior dMivamoto et al.l 
1199 It iMiyamot o et al.lll993h . All of fliis is b ehavior famil- 
iar from studies of other BHC (iHoman et al.l ,200 1, and ref- 
erences therein). Unfortunately, observation of these behav- 
iors provides little strong input to estimates of the system pa- 
rameters, i.e., mass, distance, and inclination. 'Flip-flop' be- 
havior has often been ass ociated w ith luminositi es of a few 
tens of percent (Mi yamo to^etaL|[T99 1; Miyamoto et aljl993t 
iHoman et al.ll2001l) . which is higher than we have presumed 
for the fits described here. The threshold luminosity for such 
behavior, however, is not firm, and we further note that high 
accretion rate/large distance fits to the PCA spectra failed 
04. lb . The fact that we do not see a transition to a BHC hard 
state — although we may see some indications from increases 
in the fitted disk radius at low luminosity — strongly suggest 
that the mass is > 3M0 and the distance is > lOkpc. 

Despite the fact that we do not have a firm estimate of mass 
and distance for 4U 1957-1-1 1, the relativistic disk models still 
strongly statistically prefer rapid spin solutions. Essentially, 
this is because the diskbb models fit very high tempera- 
tures - as high as 1.7 keV - with fairly low normalizations. 
We saw, however, that for the Chandra observation, or for 
the RXTE observations where the power-law/Comptonization 
component strengthens, the need for spectral hardening of the 
disk, either via a color-correction factor (hj) or rapid spin, 
was greatly reduced. Phenomenologically, a coronal compo- 
nent can mimic the effects of rapid spin. This naturally leads 
to the question of whether or not one is sure that any residual 
coronal component completely vanishes at low flux. Does one 
ever observe a "bare" disk spectrum? 

It is a rather remarkable fact that the Chandra observa- 
tion, and the low flux RXTE observations, are mostly de- 
scribed via three parameters: absorption, diskbb temper- 
ature, and diskbb normalization. Even more remarkably, 
changes in the spectrum, both in terms of amplitude and over- 
all shape, are mainly driven by variations of only one param- 
eter: diskbb temperature. The relativistic disk models, on 
the other hand, must ascribe essentially this one parameter to 
a combination of four parameters: a*, disk inclination (which 

affects the appearance of relativistic features), (M/M ), and 
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hd. In principle, any of the latter three can vary amon g obser- 
vation s. (Disk inclination can vary due to warping; see lPringlg 
11996*.) Unfortunately, there is no truly unique, sharp feature 
in the high resolution spectrum that completely breaks degen- 
eracies among these parameters. The overall magnitude of the 
diskbb temperature and shape of the 4U 1957+11 spectra, 
however, seem best reproduced with the most rapid spin mod- 
els; simply increasing accretion rates or spectral hardening 
factors is insufficient to model the spectra fully. 

We do not see the essential question as being whether or 
not 4U 1957+11 is rapidly spinning. Instead, we view the 
more observationally motivated and perhaps more fundamen- 
tal question as being: why is the characteristic disk temper- 
ature of 4U 1957+11 so high? Rapid spin is one hypothe- 
sis. Another would be that there is indeed a residual, low 
temperature corona, even at low flux. The possibly high in- 
clination of this source (« 75° to be consistent with optical 
lightcurve variations, while being consistent with the lack of 
X-ray eclipses) would mean that we are viewing the disk 
through a larger scattering depth than would be usual for 
most BHC in the soft state. We also note that there exists 
very little, self-consistent theoretical modeling of low temper- 
ature Comptonizing coronae with high temperature (l-2keV) 



seed photon temperatures. Most energetic ally bal anced , 
self-consistent models (e.g.. lDove. Wilms & B egelman 1997h 
have focused on high coronal temperatures (50-200 ke V) with 
low seed photon temperatures (100-300eV). 

If one is to take the rapid spin hypothesis for the high 
diskbb temperature at face value, then these observations 
offer something of a prediction. Statistically, they do prefer 
rapid spin, and based upon observations of other BHC we 
would have expected to observe a hard state if the faintest 
observations were < 3%LEdd- If the theoretically preferred 
value of the spectral hardening factor is indeed hj = 1.7, fu- 
ture measurements should find 4U 1957+1 1 to be an w I6M0 
black hole at w 22kpc. Whether or not this turns out to be 
the case, given the nature of 4U 1957+1 1 as perhaps the sim- 
plest, cleanest example of a BHC soft state, observations to 
independently determine this system's parameters (mass and 
distance) are urgently needed. 
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APPENDIX 

In this appendix, we describe how we model the mil d pil e up that is present in the MEG and, to lesser degree, the HEG 
data. Our model follows the discussions of lDavisI (1200 Ih and lDavisI (12003 ). although we do not attempt to arrive at an ab initio 
description of pile up. Instead, we present a simple, phenomenological description where the absolute amplitude of pile up is left 
as a fit parameter (albeit one that can be determined roughly with empirical knowledge and then frozen at a fixed value, or fitted 
with a limited range of "acceptable" values). In contrast to spectral pile up solely with the CCD detector, where piled events can 
reappear in the spectrum at higher implied energies, gratings pile up (at least in P' order spectra) can be thought of as a straight 
loss of events. Whether the piled events 'migrate' to bad grades, or are read as a single event of higher energy and are thus 
removed from the 'order sorting window', it is still a simple loss term (lDavisll2(j03h . 
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The degree of loss due to pile up scales with C,-, the number of expected incident events in a given detector region per detector 
frame integration time. The detected number of events is then reduced by a factor exp(^C,), for a suitably chosen region and 
where ^ is a 'fudge factor' of order unity. The proper detector region to consider is approximately 3 pixels in the 'cross 
dispersion' direction, by a 3-5 pixel length along the dispersion direction of the gratings arm being analyzed. (Events in a single 
frame that land in adjacent pixels will always be piled, while e vents that are removed from one another by four pixels will only 
be piled if their charge clouds extend toward one another; see lDavisll2003h Given that the peak effective area of the MEG is 
approximately twice that of the HEG, and the fact that detector pixels cover twice the wavelength range in the MEG, we expect 
the pile up to be approximately 4 times larger in the MEG. 

It is important to note that one needs to consider all events in a given detector region, i.e., all spectral orders, background events, 
contributions from the wings of the 0* order point spread function, etc., whether or not these events are otherwise rejected by 
order sorting. To partly account for this necessity, the pile up model described here uses information from the 2"^ and 3"* order 
ancillary response functions (arf s). Specifically, we create a 'convolution model' where, given an input spectrum, we define 
C,(A) as the first order arf multiplied by the unpiled model spectrum (yielding counts per detector bin). We further increase 
C,(A) by multiplying the unpiled model by the 2"'^ order arf, then shifting the 2"'' order wavelengths by a factor of two, and then 
rebinning and adding this contribution to C,(A). A similar procedure is used for the 3'^'' order contribution. We then use the input 
pile up fraction fit parameter, p j, to normalize the exponential reduction of the spectrum. Specifically, we multiply the counts per 
bin predicted for the unpiled model by exp[log(l - pf)[Ci{X)/ max(C,)]]. 

There are two important subtleties to address. First, in certain regions of the spectrum, specifically those wavelength regions 
that are dithered over chip gaps or detector bad pixels, the counts are lower not due to a low intrinsic flux, but rather due 
to a fractionally reduced exposure. The FRACEXPO column found in the arf FITS files for Chandra gratings data provides 
this information; therefore, we can approximately correct for this effect. Second, for Chandra gratings spectra, a portion of 
the detector area information is contained within the response matrix files (rmf ), which are not normalized to unity. In practice, 
before fitting the pile up model, we renormalize the gratings P' order arf s using the ISIS functions f actor_rsp, get_arf , 
and put_arf . Assuming a_id is the data index of the P' order arf, and r_id is the data index of the P' order rmf, we 
proceed as follows in ISIS: 
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It is precisely because ISIS has the abiUty to read and manipulate information from all input data files, and then easily mathe- 
matically manipulate this information to create new functionality, that we are able to define a relatively simple pile up correction 
model. 

Ideally, we should incorporate additional information, such as the variation of the line spread function (LSF; i.e., the cross 
dispersion profile) along the gratings arms, the 5-10% of the events that are typically excluded by order sorting'', etc. For these 
reasons we cannot assign an absolute estimate of the pile up fraction. However, to the extent that we have correctly captured the 
relative changes of total count rate along the dispersion direction, and have chosen a suitable peak pile up fraction, we should 
have a reasonably accurate pile up correction as a function of wavelength. More sophisticated ab initio models currently are 
under development (J. Davis, priv. comm.). 

For our data, the peak of the detected MEG count rate (including contributions from 2"'' and 3'^'' order events) occurs in the 
6-8 A region, and peaks at a value of « 0.08 counts/frame/3 pixels. Thus, making a simple first order correction since this rate 
already includes pile up, and including up to 5 pixels, we expect a peak pile up fraction in the 0.09-0.16 range. We could have 
chosen a value in this range and have frozen it in the fits, but the fitted values fall within the middle of these estimates. Likewise, 
for the HEG spectra we estimate a peak pile up fraction of 0.025-0.036. Such fractions are comparable to or smaller than existing 
systematic uncertainties in the Chandra calibration. Additionally, we found a great deal of degeneracy in any attempts to fit the 
HEG value directly, thus we froze peak the pile up fraction at a value of 0.03. 

The S-lang code that we used to define a simple gratings pile up model as an ISIS user model is presented in the electronic 
version of this manuscript. 



define simple_gpile_f it (lo, hi, par, fun) 
{ 

% Peak pileup correction goes as: 

% exp (log (1-pfrac) * [counts/max (counts) ] ) 

variable pfrac = par[0]; 

% Pileup scales with model counts from *data set* indx . . . 
variable indx = typecast (par [ 1 ], Integer_Type) ; 



' LSF variations and tlie exclusion of events that fall outside of the order sorting windows is, of course, accounted for in calculation of the gratings a r f s . 
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if ( indx == or pfrac == 0. ) 
{ 

return fun; % Quick escape for no changes . . . 

} 

% ... but the arf index could be a different number, so get that 

variable arf_indx = get_data_info (indx) . arf s; 

% Get arf information 

variable arf = get_arf (arf_indx [ ] ) ; 

% In dither regions (or bad pixel areas), counts are down not 
% from lack of area, but lack of exposure. Pileup fraction 

% therefore should scale with count rate assuming full exposure. 
% Use the arf "fracexpo" column to correct for this effect 

variable fracexpo = get_arf_info (arf_indx [ ]). fracexpo; 

if ( length (fracexpo) > 1 ) 

{ 

fracexpo [where ( fracexpo ==0)] = 1.; 

} 

else if ( fracexpo==0 ) 
{ 

fracexpo = 1; 

} 

% Rebin arf to input grid, correct for fractional exposure, and 
% multiply by "fun" to get ("corrected") model counts per bin 

variable mod_cts; 

mod_cts = fun * rebin ( lo, hi, 

arf.bin_lo, arf.bin_hi, 

arf .value/ fracexpo* (arf .bin_hi-arf .bin_lo) ) 

/ (hi-lo); 
% Go from cts/bin/s -> cts/angstrom/s 
mod_cts = mod_cts/ (hi-lo) ; 

% Use 2nd and 3rd order arfs to include their contribution. 

% Will probably work best if one chooses a user grid that extends 

% from 1/3 of the minimum wavelength to the maximum, and has 

% at least 3 times the resolution of the first order grid. 

variable mod_ord; 
if(par[2] > 0) 
{ 

indx = typecast (par [2 ], Integer_Type) ; 

arf = get_arf (indx) ; 

f racexpo=get_arf_inf o (indx) . fracexpo; 

if ( length ( fracexpo) > 1 ) 

{ 

fracexpo [where ( fracexpo ==0)] = 1.; 

} 

else if ( fracexpo==0 ) 
{ 

fracexpo = 1; 

} 

mod_ord = arf . value/f racexpo* 

rebin (arf .bin_lo, arf .bin_hi, lo, hi, fun) ; 
mod_ord = rebin (lo, hi, 2*arf .bin_lo, 2*arf .bin_hi, mod_ord) / (hi-lo) ; 



16 



Nowak et al. 



mod_cts = mod_cts+mod_ord; 

} 

if(par[3] > 0) 
{ 

indx = typecast (par [3] , Integer_Type) ; 
arf = get_arf (indx) ; 

f racexpo=get_arf_inf o (indx) . fracexpo; 

if ( length ( fracexpo) > 1 ) 

{ 

fracexpo [where ( fracexpo ==0)] = 1.; 

} 

else if( fracexpo==0 ) 
{ 

fracexpo = 1; 

} 

mod_ord = arf . value/f racexpo* 

rebin (arf .bin_lo, arf .bin_hi, lo, hi, fun) ; 
inod_ord = rebin (lo, hi, 3*arf .bin_lo, 3*arf .bin_hi, mod_ord) / (hi-lo) ; 
mod_cts = inod_cts+mod_ord; 



variable max_mod = max (mod_cts ) ; 
if (inax_inod <= 0) { max_mod = 1.; } 

% Scale maximum model counts to pfrac 

mod_cts = log (1-pf rac) *mod_cts/max_mod; 

% Return function multiplied by exponential decrease 

fun = exp(mod_cts) * fun; 

return fun; 



add_slang_f unction ( "simple_gpile" , 

["pile_frac", "data_indx", "arf2_indx", "arf3_indx"] ) ; 

set_function_category ( "simple_gpile" , ISIS_FUN_OPERATOR) ; 



